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RESUMO 


Serao desert tas brevemente neste artigo 
experiencias de pesquisa em processamento digital de 
l) deteepao de bordas por meio de testes de hipotese; 
de imagens por filtros de resposta finita ao impulso; 
trapao espacial de atributos em classif icapao mul t iespectra 1 ; 4) 
tro translacional de imagens por testes seqlienciais de hipoteses 


quatro diferentes 
imagens no INPE: 

2) interpolapao 

3) metodos de e 
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EXPERIMENTS WITH DIGITAL IMAGE PROCESSING AT INPE 


ABSTRACT 

In this article four differents research experiments 
with digital image processing at INPE will be described:. 1) edge 
detection by hypothesis testing; 2) image interpolation by finite 
impulse response filters; 3) spatial feature extraction methods in 
multiopectral classification; 4) translational image registration by 
sequential tests of hypotheses. 


1 - ABORDAGEM BAYES I ANA PARA A DETECCfiO DE BORDAS EM IMAGENS. 


Em p'rocessamento digital do fmagens e frequentemente n£ 
cessario delinear os limites entre duas regioes que tem diferentes nT 
vels de cinza, os qua i s permanecem aproximadamente constantes em cada 
regiao. Esta tarefa computac ionai de segmentagao de tmagens e usualrnen. 
te chamada detecqao de bordas. 

0 problema de detecpao de bordas em imagens tem recebido 
grande atengao nos ultimos anos. De modo gerai, entretanto, os metodos 
propostos tem sido heurTsticos, usualmente baseados em tecnicas de 
diferenciapao do tipo gradiente. A abordagem aqui proposta e baseada na 
teorta de decisao estatfstica. A imagem e modeiada pela soma de sinai 
e rufdo que sao independentes , adifcivos, gaussianos e auto-regress ivos 
em duas dimensoes. A solupao otima, em termos da teoria de decisao e_s 
tatTstica, leva a urn teste que decide entre hipoteses multiplas, com 
postas e sobrepostas. Para maiores detalhes o leitor deve consuitar 
Mascarenbas e Prado (1980). Considerando quatro "pixels" da imagem, de_ 
flnem~se sate possfveis hipoteses que envolvem tais "pixels". Denotan_ 
do por s o sinai sem rufdo e o v o slnal ruldoso, a decisao otima con_ 
siste em seiecionar o valor de i(i = 1,2,. ..,7), para o qual: 

Aj(v) = j cr(s) f(v/s) ds , (l) 
regiao i 

e maximo onde o(s) denota a funeao densidade de probabi 1 idade a priori 
do sinai s, e f(v/s) represents a funpao densidade de probabi 1 idade 
condicional de v, dado s. 

Entretanto, uma redefinigao do problema e necessaria e 
Isso leva a uma formulapao que envolve decisoes binaries sobre hipote 
ses nao-sobrepostas. Se a existencia de uma borda envolve urn limiar A 
e se os custos de decisao sao denotados por Cqj , ^iq* ^gg> esta 
decisao e tomada por: 
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onde Cjj e o custo de decldlr pela hipotese 1, quando a hipotesej e 
verdadeira. 


Obtem-se, asslm o conhecldo teste bayeslano que dec_[ 
de entre duas hipoteses compostas, o qual envolve uma razlo de verossj 
mi 1 han?a. 


0 numerador e o denominador da Expressao 2 sao computes 
dos sob as hipoteses de separabi 1 idade das matrlzes de covariancia dos 
processos que descrevem o sinal e o rufdo, bem como sob a hipotese de 
que ambos os processos sao gauss ianos. 


Para que a solu^ao do problema seja vi avel , devem-se usar 
tabelas que evitem a repeti^ao dos calculos numericos das Integrals. 
Alem disto, deve-se fazer uma nova aproximasao de modo que, ao inves de 
examinar os quatro ’'pixels" ruidosos ao tomar a decisao, apenas duas 
variaveis aleatorias sejam observadas, Assim, ao decidir borda do tipo 
1 (borda a **5°) versus nao-borda do tipo 1, ao inves de observar 
v(i,j) , v(i,j+l), v(i+l,j) e v(i+l,j + l), serao observados apehas 
v ( i , j ) e v(i,j + l) + v ( i *t* 1 , j ) + v(i+1,j-fl) v 1 onde v(i,j) represen 

3 ; ' 

ta o valor do "pixel" ruidoso na posifao ( 1 , j ) da imagem. 


0 numerador da expressao adota a forma: 


'III ds ) d S2 ds 3 ds Il 
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£ possTvel mostrar que no caso gaussiano e nas condigoes 
de alta ou baixa relagao slnal-rufdo a Expressao 3 pode ser obtida 
zendo uma Integragao dupla ao inves de uma Integragao quadrupla, isto 
e r a Expressao 3 pode ser dada por: 
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ds* f,(vj ,v'/s 1 ,s ') q | (s , ,s 1 ) , 
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que pode tambem ser especial izada no caso gaussiano, 

Afnda sob esta hipotese, tanto o numerador como o denomi 
nador tem a forma geral: 
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dx exp (~x 2 ) 
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. exp (-y 2 /2) dy. 
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A integragao interna: 


dy exp(-y 2 /2) 


( 6 ) 
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pode ser efetuada armazenando uma tabela de distribuigao gausslana, eji 
quanto a integragao de -® a +«° pode ser colocada na forma: 
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dx exp(-x 2 ) f (x) . 


(7) 


Esta integral pode ser calcula.da nume'r icamente pela 
aplicagao da formula de Gauss-Hermi te: 
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f dx oxp(~x 2 ) dx s £ H, f(x,), 
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zero do polfnomio de Hermlte do grau m, e os pesos 
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sao dados por: ■ 
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sendo os valores de x^ e dados por tabelas. 

Resultados exper Imenta I s foram obtidos com imagens de 
historia em quadrtnhos (Cebolinha), com 128 x 123 “pixels' 1 e diferen. 
tes reiapoes sinal-rufdo, e com Imagens reals de satelite, o que evj_ 
dencla um bom desempenho do algoritmo. 

2 - INTERPOLAQfiO DE IMAGENS POR FILTROS DE RESPOSTA F1NITA AO 1MPULS0. 


0 problema de interpolapio em processamento digital de 
imagens ocorre em aplicapoes como mudanpa de escala da imagem, co_r 
repao geometrica, reamostragem, etc. A abordagem aqui apresentada e 
baseada em metodos de processamento digital de sinais que representam 
uma extensao para duas dimensoes de trabalhos anteriores em uma d_j_ 
mensao, segundo um esquema separavel nas direpoes vertical e horizon_ 
tal. 

Uma classe de Interpoladores de imagem e analisada: es_ 
tas funpoes sao obtidas atraves de tecnicas de projeto de filtros de 
resposta finita ao impulse (em Ingles F.I.R.). Seu comportamento e 
comparado com os interpoladores usuais em processamento digital de 
imagens, como o do vizinho mais proximo, o bilinear e a convolupao 
-cubica'. Para maiores detaihes o leitor podera consultar Camara Neto 
e Mascarenhas (1983) e Camara Neto (1983). 


A abordagem de proccssamento digital de slnais mostra 
que o proccsso do interpolagao, que e essenclalmente um problema do v£ 
riagao da taxa de amostragem do slnal analogico original, pode ser fo_r 
mulado como um processo de filtragem passa-baixas . A restrigao de fa 
se linear sugere a implementagao por rrelo de flltros de resposta f i n_i_ 
ta, ao inves dos flltros de resposta infinlta. 

0 interpolador bidimcns ional e obtido em dots passos: um 
filtro unidlmenslonai e projetado e a interpolagao bidimensional e 
feita pela convolugao da resposta do filtro nas d! redoes horizontal e 
vertical, Tres metodos de projeto foram usados: 

a) projeto de janelamento, 

b) projeto de filtros com fase linear e ondulagao constante (em 
Ingles "equi ripple") , 

c) projeto de interpoladores com mfnimo erro medio quadrattco no 
domfnio do tempo. 

Os filtros obtidos utilizam regioes de 4 x A "pixels" para Interpolar 
um "pixel" na imagem resultante. 

No projeto por meio de janelas, uma seqUencia finita 
w(k) e utilizada para modificar o interpolador ideal de Nyquist h(k) 
(funglo sine) de modo a obter o filtro finito: 

h(k) = h(k) .w(k) . 

Entre as janelas utilizadas estao as janelas de Hamming, 
a do co-seno, a de Kaiser e a de Papoulis. Uma janela parabolica pr£ 
posta por Shlien foi tambern empregada, tendo a forma: 

,/. % . f 2k } z . N-1 N-1 

w ( k ) = i-[urrj 2~’ — \~F 


para um filtro de tamanho N. 
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Uma tccnlca mals soflsticada de projeto cnvolve a otimj_ 
zagao segundo erf ter lo dc Chebyschev, o qual conslste na minimizagao 
' do erro maximo, num intervalo de frequencia desejado. 

Flltros de Interpolagao tambem foram propostos mlnimj^ 
zando o erro medio quadratico entre as amostras Interpoladas e a sa^ 
da do interpolador Ideal(Oetken et alii (1975)).. Os parametros do fl ltro 
dependem da fungao de autocorrelagao da entrada. Foi subsequentemente 
mostrado que este procedimento de projeto poderia tambem minimizar o 
maximo erro normal Izado, se o sinal de entrada e supos^o ser llmitado 
em falxa e espectramonte piano,. Duas. medi das de. desempenho foram . _ut; i 
llzadas para avallar o erro de i nterpolagao: a perda de resol ugao e o 
erro de Interpolagao. Foi observado que as janelas simples de Shllen 
(Janela parabol lea) e a Janela do co-seno tlveram o melhor comportamento, 
mesmo se comparadas com flltros com propriedades de otimalidade. 

Para verlficar os resultados aclma, urn experlmento de 
ampliagao de escala foi feito: uma imagem teste foi reduzlda 8 vezes 
de seu tamanho original 512 x 512' para uma imagem de 6^ x 6 1 ). A im£ 
gem reduzlda foi interpolada ao seu tamanho original pelos varios m 
• terpoladores. Os parametros caiculados foram o erro e o valor absoljj 
to do erro. Os resultados obtidos corroboram as conclusSes anterlores 
de que os interpoladores projetados sob condigoes de otimalidade nao 
tern urn desempenho melhor do que aqueles projetados pelo metodo de Jja 
nelamento. Em termos visuais, as imagens geradas pelos Interpoladores 
de Shllen, sine, co-seno e convolugao cubica sao superiores aos de 
mals. Os resultados visuais podem talvez ser interpretados em termos 
da complexidade da visao humana, que certamente nao se baseia no er 
ro quadratico como criterio de fidelidade visual. 

Teen leas heuristicas tambem foram utilizadas baseadas 
em medidas da densidade local de bordas, dados que indicam o interpo 
lador mais adequado para aquela regiao. Areas que nao tern uma grande 
densidade de bordas podem ser interpoladas por urn processo mais sim 
pies. Desta maneira, uma grande economia no processo de interpolagao 
pode ser obtida. 


0 chamado operador de maxima varlagao foi utilizado para 
medir a densidade de bordas, Para uma regiao de 2 x 2 o operador de ma 
xima variagao e deflnido por: 

MV max{(|a-b| + |c-dj), ( | a-c | + |b-d|)}. 


a b 

c d 


Fig. 1 - Regiao da imagem 2x2. 

Um procedi men to heurfstlco foi desenvolvido de modo que. 
o operador MV e comparado a um llmiar. Se MV e menor do que um dado 1_J. 
miar, a interpolagao e feita pela fungao bilinear. Caso contrario, um 
interpolador de ordem mais alta como o de Schlien ou convolugao cubica 

it 

e obtido. 


Um processo adaptativo de interpolagao tambem foi utilj. 
zado no problema de mudanga de escala. Fixando a porcentagem da inte£ 
polagao bilinear em (70 ± 3)& as imagens obtidas tfnham a mesma qualj. 
dade que as anteriores, mas o tempo de computagao foi reduzido de *i0%. 

. G algoritmo adaptativo foi aplicado ao problema de re£ 
mostrar a imagem LANDSAT, que e a mostrada numa grade retangular de 
57m x 79m. Por interpolagao esta grade foi alterada para uma mais con. 
veniente, por exemplo, 50 m x 50 m. 0 bom desempenho do algoritmo foi 
confirmado neste experimento, uma vez que uma economia de 50% no tempo 
de processamento foi observada, sem nenhuma diferenga visual. 


3 - ALGUNS EXPERIMENTS COH METODOS OE EXTRAQfiO ESPACIAL OE ATRIBUTOS 
EH CLASS I FI CAC--^ HULT 1 ESPECTRAL . 

Este paragrafo descreve o uso dc novas teen leas de cxtr£ 
gao de atrlbutos espacials locals em um amblontc de classif Icogao su. 
pervlslonada. Em adlgao, a alta dimensional I dade do vetor de atribu. 
tos aumentado e contornada por um processo de selegao de atrlbutos, ba. 
seado nos dlstancias entre as dlstrlbuigoes das classes a serem discrj_ 
mlnadas. 

0 problema de classiflcagao de Imogens mul tlespectra is 
tern side abordado tradlclonalmente atraves de atrlbutos espectrais d£ 
rivados de cada canal, Todavla, a Inclusao de atrlbutos espacials pode 

t 

ser .util. 


No caso de flltragem local, a i'magem de safda e obtida 
pela convolugao da imagem original com a resposta impulsiva do fij^ 
tro especlficado pc uma mascara, isto e, uma matriz de pesos (por 
exemplo 3 x 3 ou 5 x 5). A FIgura 2 especiflca possfveis mascaras para 
a Implementagao de um filtro espaclal passa-ba ixas. A FIgura 3 apria 
senta uma aproxtmagao do operador laplaclano digital. 


(a) 



11111 

11111 

11111 

11111 

11111 


(b) 


Fig. 2 - Flltragens passa-ba ixas. 
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Fig, 3 - Aproxlmagoes digitals do operador laplaclano. 

Lima classe de flltros pao-1 Ineares pode ser Implementada 
por dols estaglos: em prlmelro lugar, uma convolugao linear e efetu£ 
da, segulda por uma fungao nao-1 Inear pontual como, por exemplo, o va 

lo 1 ' absoluto ou a fungao raiz quadrada. Claramente, flltros nao-1 Inea 

‘\‘i pddem ser deflnldos de dlferentes maneiras. Neste artlgo, um fi! 
tro nao-1 Inear denotado por vartagao fol utiilzado. Conslderando uma 
vlztnhanga 3 x 3 ao redor do "pixel" da forma: 

a b c 
d x e 
f v g h, 

a varlagao total (T.V.) e a soma da varlagao vertical (V.V.) e a varla_ 
gao horizontal (H.V.), Isto e: 

V.V. = [ a-d | + |b-x| + | c-e | + |d-f| + |x-g| + | e-h J 

H.V. ts Ja-bj + | d-x j + |f-g| + | fa — c j + |x-e| + | g-h | 

T.V. = V.V. + H.V. 



10 


No prcsente trabalho, a tarcfa de class! ficacao fo! efc 
tuada pelo critcrio tie maxima verossimi Ihanga sob a hipotese gaussiana 
(para ma lores detalhes vcr Dutra e Mascarenhas, 198^) . 

Embora util para a dlscriminagao das classes, os proce_s 
sos de cxtragao de atrlbutos espacials podem aumentar a dimens ional ld£ 
de do algoritmo de classl f Icagao, o que pode reduzir a eficiencia com 
putaclonal e exigir um numero excessivo de amostras para treinamento. 
Portanto, torna-se necessario um processo de selepao de atributos, ou 
seja, a redugao dc< dimensional idade do problema. 

Neste trabalho a distancia de Jeff reys-Matsus I ta fo! uti_ 
Uzada, sendo dada no caso gaussiano por: 


pee, 


d \lH " 2 <’ - »>• 


B = i (u -V > T (- 1 + - 2 ) (V 
8 _i _2 : — _ 


U ) + — in 
~ 2 2 


rk E + Z 

t. '— J . — 2 


l M 


E 

— 2 


V: 


onde y. e B., 1=1,2 sao o vetor media e a matriz de covartancia da 

—i —i ’ 

c 1 asse i, respect ivamente. Foram usados como criterios para selegao de 
atributos maxima distancia media e a maxima distincia minima sobre 
todos os pares' de classes. 

A extragao de atributos espaciais utilizou filtrtis lineja 
res ou nao~1 ineares. No. ultimo caso dois processos adicionais podem 
ser empregados em cascata, de uma maneira heurlstica: 


1) Um reescalonamcnto da safda do flltro nao-l Incar como uma ope 
ra$ao para melhorar o carater gausslano no qual o estagio de 

class! ficapao esta bascado. 

* 

2) . Uma filtragem linear do t!po media com duas final [dados: a) me_ 

Ihorar a,h i potese gaussiana pela combinacao linear de varl£ 
vcls aleatorias, embora algo dependentes; b) reduzlr as compo 
• nentes de alta frequencia que contem njfdo se o flltro nao-1 I 
near for do tipo derivada. 

Os filtros lineares usados.neste trabalho foram flltros 
passa-altas ou passa-ba ixas . Os filtros passa-altas extraem informapao 
sobre a Intensidade da rugosidade local e, eventualmente , sobre sua dj_ 
repao. Os f i 1 tros passa-balxas tendem a dlminuir a variancta das cla£ 
ses forgando os vetores de atributo a estar proximos dos vetores media 
de cada classe, sob a hlpotese de que os vetores pertencentes a mesma 
classe tendem a estar proximos espectral e espac ialmente. 

Considerando principalmente as limltapoes f'sicas do 
meio computacional , fo? tomada a decisao de selecionar quatro atribu_ 
tos. 


Os experimentos foram feitos com uma imagem de satelite 
Landsat-3 que cobre a area de Ribeirao Preto. Um total de quatro expei 
rimentos foram feitos. Nos tres prime! ros foram utilizados um total de 
12 atributos, de acordo com a seguinte distribuigao: 

1) Os atributos de 1 a if correspondent aos canais originais do 
Landsat. 

2) Os atributos de 5 a 8 foram obtidos pela convolugao dos canais > 
de if a 7 do Landsat com a mascara da Figura if. 


- 12 - 


3 ) Os otributos de 9 a 12 dao informagao sobre a varlagao local 
dos canals originals. Tals atrlbutos foram obcldos por difereni 
tcs processos, dependendo do conjunto de exper Imentos. 


21 

Fig. 4 - Mascara convolucional . 

Seis classes foram definidas: l) cana-de-agucar ; 2) cana. 
de-agucar nova; 3) pasto; A) agua; 5') area urbana, 6) floresta. 

No primeiro experimento, os atrlbutos de 9 a 12 foram de 
rivados dos canais originals utillzando o operador varlagao total e 
filtragem passa-baixas atraves da mascara da Figura 4 , 

No segundo expe. imento, os atrlbutos de 9 a 12 foram obi 
tldos dos canals 4 a 7 do Landsat, seguido por uma operagao de reesca. 
lonamento e filtragem passa-baixas atraves da mascara da Figura 4 . 0 
resultado da operagao de reescal onarnen to foi proporcional a raiz qua^ 
drada da entrada, Foi observado empi ricamente que tal operagao tendia 
a melhorar o carater gaussiano dos hlstogramas. 

No tercei ro experimento, os atributos de 9 a 12 foram ob^ 
tidos dos canais originals do Landsat por urn operador laplaciano rees 
calonado e suavizado. 

Foi possFvel observar' no processo de selegao de atribjj 
tos dos primeiros tres experimentos que os canais espectrais nao foram 
selecionados em qualquer ocasiao. Este fato sugere urn quarto experimen 
to, no qual os canais espectrais sao substitu nelos canais espa 
clais, de acordo com a seguinte distribuigao: 


0 1110 
11111 
11111 
11111 


suavl 
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1) Os atributos do 1 a 3 sao os canals de 4 a 6 do Landsat 
zados pcla mascara da Flgura 5. 

2) Os atributos de 4 a 6 sao os canals de 4 a 6 do Landsat suavl 
. zados pela mascara da Flgura 4. 

3) Os atributos de 7 a 9 sao os canals A, 5 e 7 do Landsat que 
passam pela operapao de variapao total, reescalonados e suavl 
zados pela mascara da Figura 4. 


4) Os atributos de 10 a 12 sao os canals 4. 5 e 6 do Landsat que 
passam pela operapao valor absolute do laplaciano, reescalona 
do e suavlzado pela mascara da Figura 4. 
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1 1 1 
1 1 1 
1 1 1 


Fig. 5 - Filtragem por media 3x3 


Os resultados obtidos mostram alguma melhoria no desemp£ 
nho medio* da classi f icapao e uma sensfvel redupao em sua confusao me 
dia pelo uso de atributos espaciais. 


Todavia, em geral, a rejeipao media (pelo uso de urn 
miar nas funpoes d i scriminantes) apresentou uma tendenc ia ao crescimeji 
to, o que pode ser considerado como nao tao danoso como cl ass i f i capao 
erronea. A area classificada se tornou mats homogp";ia e surgtu frequen_ 
temente uma regiao de rejeipao nas bordas entre classes distintas. Isso 
e dev i do ao fa to de que as medidas de rugosidade dao valores muito aj_ 
tos nessas regioes de transipao, 0 metodo proposto permite portanto a 


*Desempenho medio e definido como a peroentagem media de alas si f icapao 
oorreta para oada area (de treinamento ou teste), ponderada pelo nume 
ro de ponto de area . Definipoes analogas sao validas para a confasao 
media e rejeiodo media , sub stituindo o termo olassificaoao oorreta por 
erro de olassificaoao e de rejeiodo , respeotivamente. 


detecpao de areas com desempenho de classlf icapao usualmente pobre e 
permite rejeita-las ao inves de class! fica-las erroneamente. 

Fol observado urn aumento na dimensional idade dos dados 
que passam de dois para tres atraves da apl icapao da transformapao de 
Karhunen-Loeve. Este aumento da dimensional Idade pode ser atribuldo ao 
uso da informapao espacial. 

1* - REGISTRO TRANSLAC10NAL POR TESTES SEQtiENCIAIS DE HIPflTESES . 

0 problema de registro translaciqnal consiste em enco_n 
trar numa imagem de referenda o melhor casamento possfvel com outra 
imagem. Uma das possTveis apUcapoes destas tecnicas consiste na cla_s 
slf icapao de imagens que utilizam dados mul t i tempora i s . 

0 problema de registro translacional tern sido abordado 
nos u 1 times anos atraves de dois tipos de tecnicas, de modo geralr as 
tecnicas correlacionai s e os algoritmos seqijenciais nos'erros acumuU 
tlvos. 0 trabalho aqui descrito se insere na segunda classe de meto 
dos. A abordagem utilizada e baseada na teoria de testes sequenciais 
de -hipoteses . 


Numa primeira fase do trabalho (Mascarenhas and Pereira, 
1983), dois modelos diferentes foram utilizados: o primeiro e baseado 
na hipotese gaussiana e usa 0 fato de a variancia do erro entre duas 
Imagens tender a ser baixa no ponto de registro. 0 segundo metodo usa 
Imagens binarias derivadas das imagens originais. 0 modelo estatTstico 
para o erro acumulado resultante e uma distribuipao binominal, e a po^ 
sipao de registro e caracteri zada por uma baixa probabi 1 idade de o 
CTO binario ser unitario, Em ambos os metodos duas seqiiencias de 1J_* 
miares sao utilizadas: uma que leva a rejeipao do ponto e outra a sua 
eventual aceitapao. El as incluem o registro de uma imagem Landsat com 
suas versoes ruidosas, casamento de diferentes canais da mesma imagem 
mul tiespectral , bem. como segmentos de duas imagens obtidas em difere£ 
tes datas. 0 registro correto foi obtido, ma maioria das vezes, mesmo 
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em condipoes de balxa relapao sinal-rufdo, com um esforpo computacio 
nal modesto. . , 

Numa segunda fase (Mascarenhas and Erthal, 1984), fol 
testado com ext to o modelo exponencial para o modulo do erro entre as 
duas Imagens a s.erem registradas. Alem disso, fol utillzado um esqujj 
ma de truncamento do't'este segundo o modelo proposto por Chien e Fu 
(1966). Flnalmente, foram estabelecidas relapoes entre os modelos bino 
mlnal e gaussrano. Obleve-se a curva da probabi 1 Idade do erro binario 
no ponto de registro como funpao da relapao sinal-rufdo. 

Supondo um modelo markoviano e separavel para a funpao 
de correlapao cruzada entre os sinais das duas imagens, foram obtidas 
as curvas do erro binario como funplo do deslocamento da posipao de 
registro para varios valores dos coeficientes de correlapao nas dire_ 
poes horizontal e vertical e para varios valores da relapao sinal-ru_f 
do. Finalmente, para o modelo gaussiano, foram obtidas as curvas da 
variancia do erro em funpao do deslocamento da posiplo de registro p£ 
ra diferentes valores da relapao sinal-rufdo. 
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